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What is claimed is: 

1 . A method performed by a computer for generating a 

geometric pattern from an image having a plurality of ridges and 

_ mesh points, the method comprising the steps of: 
b 

forming a partial differential equation by transferring 
values for positions in the image to corresponding coefficients 
of the partial differential equation; 

determining simultaneous difference equations corresponding 
to the partial differential equation and the image mesh points; 
solving the simultaneous difference equations; and 
mapping the solutions of the simultaneous difference 
□ equations to respective positions on the image to determine 
y features of the image. 

Ul 2. The method of claim 1, wherein the image is a 

sS fingerprint . 

The method of claim 1, wherein the image is a facial 

The method of claim 1, wherein the image is a hand-palm 

The method of claim 1, wherein the image is an eye iris 

The method of claim 1, wherein the image is a texture 

The method of claim 1, wherein the image is an eye 
retina image. 

8. The method of claim 1, wherein the step of forming a 
partial differential equation comprises the steps of: 



3. 

20 image - 
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calculating a plurality of intrinsic properties of the image 
according to image ridge pattern; 

mapping the plurality of intrinsic properties into 
coefficients of the partial differential equation; and 

determining a boundary condition for the partial 
differential equation from the image to establish a relationship 
between properties of the image and the partial differential 
equation. 

9. The method of claim 1, wherein the step of determining 
simultaneous difference equations comprises the step of replacing 
each derivative and variable items of the partial differential 
equation with respect to each mesh point in the image. 

10. The method of claim 1, wherein the step of forming a 
partial differential equation comprises the steps of: 

determining initial conditions for the partial differential 
equation; and 

determining boundary condition for the partial differential 
equation. 

11. The method of claim 10, wherein the step of determining 
initial conditions comprises the steps of: 

normalizing the image to reduce variations in gray-level 

values along ridges of the image ; 

estimating property values of the image; and 

mapping the estimated property values into weight 

coefficients of the partial differential equation. 

12. The method of claim 10, wherein the step of determining 
boundary condition comprises the steps of: 

drawing a close boundary within the image; and 
setting boundary condition on the drawn boundary. 
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13. The method of claim 11, wherein the step of normalizing 
comprises the steps of: 

determining M, mean of the gray-level in a region R by 
computing 

M=(1/N)SX F(IJ) 
(U)eR 

where N is total number pixels in the region R; F(I, J) is 
gray value of the image at point (I, J) ; 
determining V, variance of the region R by computing 

V=(l/N) 2 2 (F(I,J)-M(F))*(F(I,J)-M(F)) 

(2b) ; and 

(U)eR 

determining a normalized region R by computing 

R(IJ)=m+sqrt((v*(F(IJ)-M)*(F(IJ)-M))/V),if (U)>M; 
R(I J)=m-sqrt((v*(F(IJ)-M)*(F(U)-M))/V), otherwise 

where m and v are the desired mean and variance values, 
respectively. 

14. The method of claim 11, wherein the step of estimating 
property values comprises the steps of: 

dividing a region R into blocks of size b*b as B(k); 
computing gradients at each pixel in R as 

ax(IJ)=(pl*F(I-d 9 J)+p2*F(IJ)+p3*(F(I+dJ))/p, 
a y (IJ)=(pl*F(IJ-d)+ P 2*F(U)+p3*(F(IJ+d))/p 

where dx (I, J) and 9y (I, J) are the gradient magnitude 
in x and y directions at pixel (I, J) of the image 
respectively, pi, p2 are positive numbers, p2 is negative 
number, and d is a constant expressed as step of the 
gradients p = pl+p2+p3; 
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estimating local orientation of each block (B(k)) centered 
at pixel (I, J) by computing 

£x(U)= S 1 2*ax(u,v)*3y(u,v), ( ^ 
(u,v)GB(k) 

Cy(U)= X S 03x(u,v)*3ay(u,v)), ( 3c ) 
(u,v)eB(k) 

0(I,J)=(l/2)atan{|x(I,J)/Cy(U)} ( 3d > 

Where 0(I,J) is an estimate of the local ridge 
orientation at the block centered at pixel (I, J) ; 

computing ridge oriental vector as 

P=(l/n)2cos(2*0(I,J)) 

(4a) 

(U)eR 

Q=(l/n)£sin(2*0(I,J)) ^ 
(U)eR 

where P and Q are first and second components of the 
ridge oriental vector, respectively, and n is the 
total number of pixels calculated at the local region 
R; 

for each block centered at pixel (I, J), computing the 
minimal value and maximal value the pixel; 

for each block centered at pixel (I, J), computing a 
sequence of pixels seq(I, J) that take minimal and maximal value 
along the direction (a, b) , where (a, b) is orthogonal vector of 
the main oriental vector (p, q) ; 
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computing freq(I, J), frequency of seq(I, J) at each block 
centered at pixel (I, J) according to the differential value 
between connected elements in seq(I, J) ; and 

estimating a local ridge frequency W by computing 

W=(l/K)SIfreq(u,v) 

. (5) 

(u,v)G wnd(I,J) 

15. The method of claim 11, wherein the step of mapping 
comprises the step of determining weight coefficients, Al, A2 , 
A3, A4, A5 and A6 of a partial differential equation 

A1 ^ +A2 ^ + A3^ + A4^ + A 5 H +A6 * U=0 
3X 2 3X3Y 3Y 2 3X 3Y 

by computing 

A1=P*P*(P*P+Q*Q)*W*W, (6a) 

A2=2*(sqrt(u*u-P*P*W*W)*sqrt(v-Q*Q*W*W))/W, ( 6b) 

A3=Q*Q(P*P+Q*Q)*W*W, ( 6c ) 

A4=u*q+v, (6d) 

A5=-V*p-U, and (6e) 

A6=a*(P*P+Q*Q)+b ( 6 f ) 

where a, b, u, v are constants. 



16. The method of claim 12, wherein the step of setting 
boundary condition comprises the steps of: 

denoting B as a discrete boundary in a region R; and 
computing 



U 



=Bl(x,y), 
(x,y)eS 



au 



as 



=B2(x,y) 
(x,y)eS 



(7) 



where S is a continuous boundary defined on the discrete 
boundary B, and Bl (x, y) and B2 (x, y) are the continuous 
function and dif f erentiable function defined on the 
boundary S, respectively. 
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17. The method of claim 1, wherein the step of determining 
simultaneous difference equations comprises the steps of: 

integralizing the image to produce a group of integral 
points within a region R and an integral boundary IB; 

discretizing the image based on the mesh points for 
numerating the partial differential equation and a boundary 
condition; and 

transforming the discretizied image by solving each mesh 
point in the image to determine the simultaneous difference 
equations . 

18. The method of claim 17, wherein the step of 
integralizing comprises the steps of: 

denoting two directions of the coordinate axes of the 
fingerprint image as X-direction and Y-direction; and 

along the X-direction and Y-direction, computing integral 
points at a desired step length H as follow 

X(I) = XO + I*H, 1 = 0, 1, 2, , W(F), (8a) 

Y(J) = Y0 + J*H, J = 0, 1, 2, , H(F), (8b) 

where, (XO, Y0) is top left point of the image, W(F) is the 
width of the image and H(F) is the height of the image. 

19. The method of claim 17, wherein the step of 
discretizing comprises the steps of: 

computing the derivatives in the partial differential 

equation 

a 2 u a 2 u AO a 2 u A „au AC au n M . 
A1 8F +A2 a^ +A3 a^ +A4 ^ +A5 ^ +A6 U=0 

with respect to each inner mesh point as 

— «[U(X+H,Y)-U(X,Y)]/H (9a) 

ax 

9U =[U(X,Y+H)-U(X,Y)]/H (9b) 



aY 
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a 2 u 



ax 2 
a 2 u 



dY 

d 2 U 



[U(X+H,Y)-2*U(X,Y)+U(X-H,Y)]/H*H) ( 9c ) 

2 [U(X,Y+H)-2*U(X,Y)+U(X,Y-H)]/(H*H) (9d) 



=[U(X+H,Y+H)-U(X+H,Y)-U(X,Y+H)+U(X,Y)]/(H*H) (9e) 



axax 

where (X, Y) is inner mesh point in region R, IMP(R); 
discretizing the boundary condition 

Bl(x,y), 3U 



u 



(x,y) E S as 



= B2(X ' y) .-ana 



(x,y)GS 

combining the numerical derivatives of the partial 
differential equation and numerical boundary condition. 

20. The method of claim 19, wherein the step of 
discretizing the boundary condition comprises the steps of: 

replacing the continuous function U(X, Y) in the boundary 
condition with discrete function F(I, J) , wherein (I, J) is inner 
mesh point of a region in the image; 

replacing the continuous function Bl(x, y) in the boundary 
condition with a numerical function according to 

Dl(X,Y)=fl*F(X,Y)+f2, (X,Y)E IB (10) 

where fl and f2 are constants that are predetermined 

according to brightness and contrast of the image, and F(X, 

Y) is the gray value at point (X, Y) on the integral 

boundary IB; and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X 9 Y)=fl*[F(Xl,Yl)-F(X,Y)]/h, (X,Y) £ IB 

h=sqrt((Xl-X)*(Xl-X)+CYl-Y)*(Yl-Y)); 
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where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

21. The method of claim 1, wherein the step of mapping the 
solutions comprises the steps of: 

for each region R, computing minimum value min (R) and 
maximum value max (R) in the solution; 

computing ratio value r (R)= 255/ (max (R) - min (R) ) ; 

for each point (I, J) in the region R, assigning v(I, J) * 
3 as the solution at point (I, J) ; 

computing W(I, J) = r(R)*(v(I, J)-min(R)) as the mapping of 
w(I, J) into gray level byte at the position (I, J); and 

enhancing the region R is by placing the value W(I, J) at 
position (I, J) . 

22. A method performed by a computer for extracting 
features from an image, the method comprising the steps of: 

establishing a mathematical model according to regional 
conditions in the image; 

converting the mathematical model into numerical equations ; 
solving the numerical equations; and 

transferring the solutions of the numerical equations to 
respective regions of the image. 

23. The method of claim 22, wherein the step of 
establishing a mathematical model comprises the steps of: 

forming a partial differential equation; 

calculating a plurality of intrinsic properties of the image 
according to image ridge pattern; 

mapping the plurality of intrinsic properties into 
coefficients of the partial differential equation; and 

determining a boundary condition for the partial 
differential equation from the image to establish a relationship 
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between properties of the image and the partial differential 
equation. 

24. The method of claim 22, wherein the step of 
establishing a mathematical model comprises the steps of: 

forming a partial differential equation; 

determining initial conditions for the partial differential 
equation; and 

determining boundary condition for the partial differential 
equation. 

25. The method of claim 24, wherein the step of determining 
initial conditions comprises the steps of: 

normalizing the image to reduce variations in gray-level 
values along ridges of the image; 

estimating property values the image; and 

mapping the estimated property values into weight 
coefficients of the partial differential equation. 

26. The method of claim 24, wherein the step of determining 
boundary condition comprises the steps of: 

drawing a close boundary within the image; and 
setting boundary condition on the drawn boundary. 

27. The method of claim 25, wherein the step of normalizing 
comprises the steps of: 

determining M, mean of the gray-level in a region R by 
computing 

M=(1/N)2S F(U) 
(U)eR 

where N is total number pixels in the region R; F(I, J) is 
gray value of the image at point (I, J) ; 
determining V, variance of the region R by computing 
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V=(l/N) 2 £ (F(I,J)-M(F))*(F(I,J)-M(F)) 
(U)eR 



(2b) ; and 



determining a normalized region R by computing 



R(I,J)=m+sqrt((v*(F(I,J)-M)*(F(I,J)-M))A^),if (I,J)>M; 
R(I,J)=m-sqrt((v*(F(I,J)-M)*(F(I,J)-M))A^), otherwise 



(2c) 



where m and v are the desired mean and variance values, 
respectively . 

28. The method of claim 25, wherein the step of estimating 
property values comprises the steps of: 

dividing a region R into blocks of size b*b as B(k); 
computing gradients at each pixel in R as 



where dx (I, J) and 3y (I, J) are the gradient magnitude 
in x and y directions at pixel (I, J) of the image 
respectively, pi, p2 are positive numbers, p2 is negative 
number, and d is a constant expressed as step of the 
gradients p = pl+p2+p3; 

estimating local orientation of each block (B(k)) centered 
at pixel (I, J) by computing 



ax(I,J)=(pl*F(I-d,J)+p2*F(I,J)+p3*(F(I+d,J))/p, 
9y(I,J)=(pl*F(IJ-d)+p2*F(I,J)+p3*(F(I,J+d))/p 



(3a) 



£x(I,J)= X 1 2*dx(u,v)*ay(u,v), 

(u,v)eB(k) 
Cy(U)= X X 0dx(u,v)*3dy(u,v)), 

(u,v)GB(k) 
d (I, J)=( 1 /2)atan { £ x(I, J)/ £ y(IJ) } 



(3b) 



(3c) 



(3d) 
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where 9(1,1) is an estimate of the local ridge 
orientation at the block centered at pixel (I, J) ; 

computing ridge oriental vector as 
PKl/n)Xcos(2*6>(I,J)) 

(4a) 

(U)gR 

Q=(l/n)Esin(2*0(IJ)) (4b) 
(U)eR 

where P and Q are first and second components of the 
ridge oriental vector, respectively, and n is the 
total number of pixels calculated at the local region 
R; 

for each block centered at pixel (I, J), computing the 
minimal value and maximal value the pixel; 

for each block centered at pixel (I, J), computing a 
sequence of pixels seq(I, J) that take minimal and maximal value 
along the direction (a, b) , where (a, b) is orthogonal vector of 
the main oriental vector (p, q) ; 

computing freq(I, J), frequency of seq(I, J) at each block 
centered at pixel (I, J) according to the differential value 
between connected elements in seq(I, J) ; and 

estimating a local ridge frequency W by computing 

W=(l/K)£Xfreq(u,v) 

. (5) 

(u,v)e wnd(IJ) 

29. The method of claim 25, wherein the step of mapping 
comprises the step of determining weight coefficients, Al, A2 , 
A3, A4, A5 and A6 of a partial differential equation 
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by computing 
A1=P*P*(P*PK>*Q)*W*W, 

A2=2*(sqrt(u*u-P*P*W*W)*sqrt(v-Q*Q*W*W))/W, 

A3=Q*Q(P*P+Q*Q)*W*W, 

A4=u*q+v, 

A5=-v*p-U, and 

A6=a*(P*P+Q*Q)+b 



(6a) 



(6b) 
(6c) 
(6d) 
(6e) 
(6f) 



where a, b, u, v are constants. 

30. The method of claim 26, wherein the step of setting 
boundary condition comprises the steps of: 

denoting B as a discrete boundary in a region R; and 
computing 



where S is a continuous boundary defined on the discrete 
boundary B, and Bl(x, y) and B2 (x, y) are the continuous 
function and dif f erentiable function defined on the 
boundary S, respectively. 

31. The method of claim 22, wherein the step of converting 
the mathematical model comprises the steps of: 

determining simultaneous difference equations; 

integralizing the image to produce a group of integral 
points within a region R and an integral boundary IB; 

discretizing the image based on mesh points for numerating 
the partial differential equation and a boundary condition; and 

transforming the discretizied image by solving each mesh 
point in the image to determine the simultaneous difference 
equations . 

32. The method of claim 31, wherein the step of 
integralizing comprises the steps of: 



=Bl(x,y), 
(x,y)eS 



au =B2(x,y) 



(7) 



9s (x,y) g S 
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denoting two directions of the coordinate axes of the 
fingerprint image as X-direction and Y-direction; and 

along the X-direction and Y-direction, computing integral 
points at a desired step length H as follow 

X(I) = XO + I*H, 1 = 0, 1, 2, , W(F), (8a) 

Y(J) = Y0 + J*H, J = 0, 1, 2, , H (F) , (8b) 

where, (XO, Y0) is top left point of the image, W(F) is the 
width of the image and H(F) is the height of the image. 

33. The method of claim 31, wherein the step of 
discretizing comprises the steps of: 

computing the derivatives in the partial differential 
equation 

with respect to each inner mesh point as 

— =[U(X+H,Y)-U(X,Y)]/H (9a) 
9X 

— - [U(X, Y+H)-U(X,Y)]/H ( 9b ) 

dY 

« [U(X+H,Y)-2*U(X, Y)+U(X-H,Y)]/H*H) (9c) 

dX 2 

^-«[U(X,Y+H>2*U(X,Y)+U(X,Y-H)]/(H*H) ( 9d) 

dY 

-i^L- [U(X+H,Y+H)-U(X+H,Y)-U(X,Y+H)+U(X,Y)]/(H*H) (9e) 
wnere (X, Y) is inner mesh point in region R, IMP(R); 

discretizing the boundary condition 
=Bl(x,y), 9U 



U 



(x,y) e S ds 



=B2(X ' y) ;and 

(x,y)eS 
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combining the numerical derivatives of the partial 
differential equation and numerical boundary condition. 

34. The method of claim 33, wherein the step of 
discretizing the boundary condition comprises the steps of: 

replacing the continuous function U(X, Y) in the boundary 
condition with discrete function F (I, J) , wherein (I, J) is inner 
mesh point of a region in the image; 

replacing the continuous function Bl(x, y) in the boundary 
condition with a numerical function according to 

Dl(X,Y)=fl*F(X,Y)+f2, (X,Y)E IB (10) 

where fl and f2 are constants that are predetermined 

according to brightness and contrast of the image, and F(X, 

Y) is the gray value at point (X, Y) on the integral 

boundary IB; and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X,Y)=fl*[F(Xl 9 Yl)-F(X,Y)]/h ? (X,Y) £ IB 

h=sqrt((Xl-X)*(Xl-X)+(Yl-Y)*(Yl-Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

35. The method of claim 22, wherein the step of 
transferring the solutions comprises the steps of: 

for each region R in the image, computing minimum value min 
(R) and maximum value max (R) in the solution; 

computing ratio value r (R)= 255/ (max (R) - min (R) ) ; 

for each point (I, J) in the region R, assigning v(I, J) * 
3 as the solution at point (I, J); 

computing W(I, J) - r(R)*(v(I / J)-min(R)) as the mapping of 
w(I, J) into gray level byte at the position (I, J); and 
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enhancing the region R is by placing the value W(I, J) at 
position (I, J) . 

36. The method of claim 22, wherein the image is one or 
more of a fingerprint image, a facial image, a hand-palm image, 
an e y e iris, an eye retina, and a texture image. 

37. A digital signal processor (DSP) having stored thereon 
a set of instructions including instructions for generating 
geometric pattern from an image having a plurality of ridges and 
mesh points, when executed, the instructions cause the DSP to 
perform the steps of: 

forming a partial differential equation by transferring 
values for positions in the image to corresponding coefficients 
of the partial differential equation; 

determining simultaneous difference equations corresponding 
to the partial differential equation and the image mesh points; 

solving the simultaneous difference equations; and 

mapping the solutions of the simultaneous difference 
equations to respective positions on the image to determine 
features of the image. 

38. The DSP of claim 37, wherein the step of forming a 
partial differential equation comprises the steps of: 

determining initial conditions for the partial differential 
equation; and 

determining boundary condition for the partial differential 
equation. 

39. The DSP of claim 38, wherein the step of determining 
initial conditions comprises the steps of: 

normalizing the image to reduce variations in gray-level 
values along ridges of the image; 

estimating property values the image; and 
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mapping the estimated property values into weight 
coefficients of the partial differential equation. 

40. The DSP of claim 38, wherein the step of determining 
boundary condition comprises the steps of: 

drawing a close boundary within the image; and 
setting boundary condition on the drawn boundary. 

41. The DSP of claim 39, wherein the step of normalizing 
comprises the steps of: 

determining M, mean of the gray-level in a region R by 
computing 

M=(1/N)SS F(IJ) 

(2a) 

CU)eR 



where N is total number pixels in the region R; F(I, J) is 
gray value of the image at point (I, J) ; 
determining V, variance of the region R by computing 

V=(l/N) 1 1 (F(IJ)-M(F))*(F(IJ)-M(F)) 

(2b) ; and 

(U)€R 

determining a normalized region R by computing 
R(U)=m+sq^ 

R(I ,J)=m-sqrt((v*(F(I ,J)-M)*(F(I,J)-M))A^), otherwise 

where m and v are the desired mean and variance values, 
respectively. 

42. The DSP of claim 39, wherein the step of estimating 
property values comprises the steps of: 

dividing a region R into blocks of size b*b as B(k); 
computing gradients at each pixel in R as 
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ax(I,J)=(pl*F(I-d,J)+p2*F(I,J)+p3*(F(I+d,J))/p, 
a y a,J)=(pl*F(I,J-d)+p2*F(I,J)+p3*(F(U+d))/p 

where 8x (I, J) and dy (I, J) are the gradient magnitude 
in x and y directions at pixel (I, J) of the image 
respectively, pi, p2 are positive numbers, p2 is negative 
number, and d is a constant expressed as step of the 
gradients p = pl+p2+p3 ; 

estimating local orientation of each block (B(k)) centered 
at pixel (I, J) by computing 

£x(IJ)= £ S 2*dx(u,v)*dy(u,v), 

(3b) 

(u,v)GB(k) 

£y(U)= S X 03x(u,v)*3ay(u,v)), { 3 c } 

(u,v)eB(k) 

0 (I, J)=( 1 /2)atan { | x(I, J)/ £ y(U) } (3d) 

Where 0(I,J) is an estimate of the local ridge 
orientation at the block centered at pixel (I, J) ; 

computing ridge oriental vector as 

P=(l/n)2cos(2*0(U)) (4a) 
(U)gR 

Q=(l/n)i:sin(2*0(U)) (4b) 
(U)eR 

Where P and Q are first and second components of the 
ridge oriental vector, respectively, and n is the 
total number of pixels calculated at the local region 
R; 

for each block centered at pixel (I, J), computing the 
minimal value and maximal value the pixel ; 
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for each block centered at pixel (I, J), computing a 
sequence of pixels seq(I, J) that take minimal and maximal value 
along the direction (a, b) , where (a, b) is orthogonal vector of 
the main oriental vector (p, q) ; 

computing freq(I, J), frequency of seq(I, J) at each block 
centered at pixel (I, J) according to the differential value 
between connected elements in seq(I, J) ; and 

estimating a local ridge frequency W by computing 
W=(l/K)SSfreq(u,v) 

(5) 

(u,v)e wnd(I,J) 

43. The DSP of claim 39, wherein the step of mapping 
comprises the step of determining weight coefficients, Al, A2 , 
A3, A4, A5 and A6 of a partial differential equation 

A1 ^ +A2 ^ +A3 ^ + A4™ + A5#Wu=0 
ax 2 3Y 2 8X 3Y 

by computing 



A1=P*P*(P*P+Q*Q)*W*W, (6a) 

A2=2*(sqrt(u*u-P*P*W*W)^^^ (6b) 

A3=Q*Q(P*P+Q*Q)*W*W, ( 6c ) 

A4=u*q+v, (6d) 

A5=-V*p-U ? and (6e) 

A6=a*(P*P+Q*Q)+b (6f ) 
where a, b, u, v are constants. 



44. The DSP of claim 40, wherein the step of setting 
boundary condition comprises the steps of: 

denoting B as a discrete boundary in a region R; and 
computing 
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u 



=Bl(x,y), 
(x,y)eS 



au =B2(x,y) 
3s (x,y)eS 



(7) 



where S is a continuous boundary defined on the discrete 
boundary B, and Bl (x, y) and B2 (x, y) are the continuous 
function and dif f erentiable function defined on the 
boundary S, respectively. 

45. The DSP of claim 37, wherein the step of determining 
simultaneous difference equations comprises the steps of: 

integralizing the image to produce a group of integral 
points within a region R and an integral boundary IB; 

discretizing the image based on the mesh points for 
numerating the partial differential equation and a boundary 
condition; and 

transforming the discretizied image by solving each mesh 
point in the image to determine the simultaneous difference 
equations . 

46. The DSP of claim 45, wherein the step of integralizing 
comprises the steps of: 

denoting two directions of the coordinate axes of the 
fingerprint image as X-direction and Y-direction; and 

along the X-direction and Y-direction, computing integral 
points at a desired step length H as follow 



where, (XO, YO) is top left point of the image, W(F) is the 
width of the image and H(F) is the height of the image. 

47. The DSP of claim 45, wherein the step of discretizing 
comprises the steps of: 

computing the derivatives in the partial differential 
equation 



X(I) = XO + I*H, 
Y(J) = YO + J*H, 



1 = 0, 1, 2, 
J = 0, 1, 2, 



W(F) , 
H(F) , 



(8a) 
(8b) 
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A1 ^ +A2 i!lL +A 3^ + A4|i + A5^ + A6.U=0 
3X 2 3X3Y 8Y 2 3X SY 

with respect to each inner mesh point as 

— «[U(X+H,Y)-U(X,Y)]/H ( 9a) 

r)TT 

— « [U(X,Y+H)-U(X, Y)]/H ( 9b ) 

a 2 u 



ax 2 
a 2 u 



9Y 2 

a 2 u 



- [U(X+H,Y)-2*U(X,Y)+U(X-H, Y)]/H*H) (9c) 
[U(X,Y+H)-2*U(X,Y)+U(X,Y-H)]/(H*H) ( 9d) 



: [U(X+H,Y+H)-U(X+H,Y)-U(X,Y+H)+U(X,Y)]/(H*H) (9e) 



wiiere (X, Y) is inner mesh point in region R, IMP(R) 
discretizing the boundary condition 

=Bl(x,y), au 



u 



(x,y)sS as 



=B2(X ^ ;and (7) 



(x,y)eS 

combining the numerical derivatives of the partial 
differential equation and numerical boundary condition. 

48. The DSP of claim 47, wherein the step of discretizing 
the boundary condition comprises the steps of: 

replacing the continuous function U(X, Y) in the boundary 
condition with discrete function F (I, J), wherein (I, J) is inner 
mesh point of a region in the image; 

replacing the continuous function Bl (x, y) in the boundary 
condition with a numerical function according to 

Dl(X,Y)=fl*F(X,Y)+f2, (X,Y)eIB (10) 
where fl and f2 are constants that are predetermined 
according to brightness and contrast of the image, and F(X, 
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Y) is the gray value at point (X, Y) on the integral 
boundary IB; and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X 9 Y)=fl*[F(Xl,Yl)-F(X,Y)]/h, (X,Y)eIB 
h=sqrt((Xl -X)*(X 1 -X)+(Y 1 -Y)*(Y 1 -Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

49. The DSP of claim 37, wherein the step of mapping the 
solutions comprises the steps of: 

for each region R, computing minimum value min (R) and 
maximum value max (R) in the solution; 

computing ratio value r (R)= 255/ (max (R) - min (R) ) ; 

for each point (I, J) in the region R, assigning v(I, J) * 
3 as the solution at point (I, J) ; 

computing W( I, J) = r(R)*(v(I / J)-min(R)) as the mapping of 
w(I, J) into gray level byte at the position (I, J) ; and 

enhancing the region R is by placing the value W(I, J) at 
position (I, J) . 

50. The DSP of claim 37, wherein the image is one or more 
of a finger print, a facial image, a hand-palm image, an eye 
iris, an eye retina, and a texture image. 

51. A method performed by a computer for biometric image 
processing, the DSP comprising the steps of: 

establishing numerical relationship between visual 
appearance of the biometric image; and 

approximating solutions of a partial differential equation 
with a boundary condition according to the established numerical 
relationship to determine features of the biometric image. 
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52. The method of claim 51, wherein the step of 
establishing numerical relationship comprises the steps of: 

forming a partial differential equation; 

calculating a plurality of intrinsic properties of the image 
according to image ridge pattern; 

mapping the plurality of intrinsic properties into 
coefficients of the partial differential equation; and 

determining a boundary condition for the partial 
differential equation from the image to establish a relationship 
between properties of the image and the partial differential 
equation. 

53. The method of claim 52, wherein the step of forming a 
partial differential equation comprises the steps of: 

determining initial conditions for the partial differential 
equation; and 

determining boundary condition for the partial differential 
equation. 

54. The method of claim 53, wherein the step of determining 
initial conditions comprises the steps of: 

normalizing the image to reduce variations in gray-level 
values along ridges of the image; 

estimating property values the image; and 

mapping the estimated property values into weight 
coefficients of the partial differential equation. 

55. The method of claim 53, wherein the step of determining 
boundary condition comprises the steps of: 

drawing a close boundary within the image; and 
setting boundary condition on the drawn boundary. 

56. The method of claim 51, wherein the step of 
approximating solutions comprises the steps of: 
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integralizing the image to produce a group of integral 
points within a region R and an integral boundary IB; 

discretizing the image based on mesh points for numerating 
a partial differential equation and a boundary condition; 

transforming the discretizied image by solving each mesh 
point in the image to determine the simultaneous difference 
equations; 

replacing the continuous function U(X, Y) in the boundary 
condition with discrete function F(I, J), wherein (I, J) is inner 
mesh point of a region in the image; 

replacing the continuous function Bl(x, y) in the boundary 
condition with a numerical function according to 

Dl(X,Y)=fl*F(X,Y)+f2, (X,Y)e IB (10) 
where fl and f2 are constants that are predetermined 
according to brightness and contrast of the image, and F(X, 
Y) is the gray value at point (X, Y) on the integral 
boundary IB; and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X,Y)=fl*[F(Xl ? Yl)-F(X,Y)]/h 9 (X,Y) £ IB 

h=sqrt((Xl-X)*(Xl-X)+(Yl-Y)*(Yl-Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

57. The method of claim 56, wherein the step of 
integralizing comprises the steps of: 

denoting two directions of the coordinate axes of the 
fingerprint image as X-direction and Y-direction; and 

along the X-direction and Y-direction, computing integral 
points at a desired step length H as follow 

X(I) = XO + I*H, 1 = 0, 1, 2, , W(F), (8a) 

Y(J) = Y0 + J*H, J - 0, 1, 2, , H(F), (8b) 
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where, (XO, YO) is top left point of the image, W(F) is the 
width of the image and H(F) is the height of the image. 

58. The method of claim 56, wherein the step of 
discretizing comprises the steps of: 

computing the derivatives in the partial differential 
equation 

with respect to each inner mesh point as 

au 



ax 
au 



av 

a 2 u 



=[U(X+H,Y)-U(X,Y)]/H ( 9a ) 

«[U(X,Y+H)-U(X,Y)]/H (9b) 



ax 2 
a 2 u 



=[U(X+H,Y)-2*U(X,Y)+U(X-H,Y)]/H*H) ( 9c ) 



aY 2 
a 2 u 



~[U(X,Y+H)-2*U(X,Y)+U(X,Y-H)]/(H*H) (9d) 



=[U(X+H,Y+H)-U(X+H,Y)-U(X,Y+H)+U(X,Y)]/(H*H) (9e) 



where (X, Y) is inner mesh point in region R, IMP(R); 
discretizing the boundary condition 

=Bl(x,y), au 



U 



=B2(X ^ ;and 

(x,y)eS 



(x,y) g S as 

combining the numerical derivatives of the partial 
differential equation and numerical boundary condition. 

59. The method of claim 56, wherein the step of 
discretizing the boundary condition comprises the steps of: 
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replacing the continuous function U(X, Y) in the boundary 
condition with discrete function F(I, J) , wherein (I, J) is inner 
mesh point of a region in the image; 

replacing the continuous function Bl (x, y) in the boundary 
condition with a numerical function according to 

Dl(X,Y)=fl*F(X,Y)+f2, (X,Y)e IB (10) 
where fl and f2 are constants that are predetermined 
according to brightness and contrast of the image, and F(X, 
Y) is the gray value at point (X, Y) on the integral 
boundary IB; and 

replacing the continuous function B2 (x, y) in the boundary 
condition with a numerical function according to 

D2(X > Y)=fl*[F(Xl,Yl)-F(X,Y)]/h, (X,Y)eIB 
h=sqrt((Xl-X)*(Xl-XH(Yl-Y)*(Yl-Y)); 

where (XI, Yl) is an integral point on IB selected as the 
next adjacent point along the boundary line IB. 

60. The method of claim 56, further comprising the step of 
mapping the solutions of the partial differential equation to 
respective positions on the image to determine features of the 
image . 

61. The method of claim 60, wherein the step of mapping the 
solutions comprises the steps of: 

for each region R, computing minimum value min (R) and 
maximum value max (R) in the solution; 

computing ratio value r (R) = 255/ (max (R) - min (R) ) ; 

for each point (I, J) in the region R, assigning v(I, J) * 
3 as the solution at point (I, J); 

computing W(I, J) = r(R)*(v(I / J)-min(R)) as the mapping of 
w(I, J) into gray level byte at the position (I, J); and 
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5 



10 



Hjl5 



U20 



25 



30 



enhancing the region R is by placing the value W(I, J) at 
position (I, J) . 

62. The method of claim 51, wherein the image is one or 
more of a finger print, a facial image, a hand-palm image, an eye 
iris, an eye retina, and a texture image. 
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